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We apply a new Kubo-Greenwood type formula combined with a generalized Feynman diagram- 
matic technique to report a first principles calculation of the thermal transport properties of dis- 
ordered Yei-xCvx alloys. The diagrammatic approach simplifies the inclusion of disorder-induced 
scattering effects on the two particle correlation functions and hence renormalizes the heat current 
operator to calculate configuration averaged lattice thermal conductivity and diffusivity. The ther- 
mal conductivity k,{T) in the present case shows an approximate quadratic T-dependence in the low 
temperature regime (T < 20 K), which subsequently rises smoothly to a T- independent saturated 
value at high T. A numerical estimate of mobility edge from the thermal diffusivity data yields 
the fraction of localized states. It is concluded that the complex disorder scattering processes, in 
force-constant dominated disorder alloys such as Fe-Cr, tend to localize the vibrational modes quite 
significantly. 

PACS numbers: 66.70.-f, 66.30.Xj,63.20.Pw 



The study of phonon excitations and the associated 
thermal transport properties is an important field of re- 
search in disordered alloys. In certain materials, dis- 
order mediated scattering can shrink the typical mean 
free path (MFP) of phonons to such a level that wave- 
length and MFP no longer remain sharp concepts, and 
the usual textbook phonon gas model for thermal con- 
ductivity breaks down. From the theoretical perspective, 
the development of a reliable quantum mechanical theory 
to predict such properties in random alloys is a difhcult 
task mainly because of two problems : (i) one needs a 
microscopic description of inter-atomic force constants 
with an intrinsic off-diagonal disorder and (ii) one has to 
configuration average a two-particle correlation function 
within a Kubo-type formula. The effects of dominant off- 
diagonal force constant disorder in alloys can be quite un- 
usual, as we have shown earlier.^ Most theories of thermal 
transport, developed in the past few decades however, are 
either based on the single-site coherent potential approx- 
imation (CPA),^ the perturbation-based approach simu- 
lating the Peicrls-Boltzmann equation (PBE)"^ or atom- 
istic models with a large unit cell and periodic boundary 
conditions.^ CPA, being a single-site mean-field approxi- 
mation, is inadequate for treating multi-site off-diagonal 
disorder arising out of force constants and is reported, for 
example, to inadequately explain experimental life-time 
data on simple Ni-Pt alloys.^ The perturbative simula- 
tion approach, although rigorously derived, is limited in 
applicability to model lattices alone and has not been 
tested on realistic materials. The atomistic models are 
computationally expensive due to the large unit-cell size, 
non-self-consistent and suffer from the finite size errors. 

In a recent paper^ we have developed a theoretical 
approach to calculate the configuration averaged lattice 



thermal conductivity and diffusivity for random alloys. 
This formalism combined a Kubo-Greenwood approach 
with a generalized Feynman diagrammatic technique to 
explicitly incorporate the effect of disorder induced scat- 
tering. We showed that disorder scattering renormalizes 
both the phonon propagators as well as the heat cur- 
rents. These corrections are related to the self-energy 
and vertex corrections. Unlike the single-site CPA, this 
approach explicitly takes into account the fluctuations 
in masses (diagonal) , force constants and heat currents 
(off-diagonal disorder) between different ion-cores and in- 
corporates the sum rule relating the diagonal element of 
the force constant to the off-diagonal ones. 

In the present paper, we combine this theoretical 
approach with a first-principles Quantum-ESPRESSO 
(QE) calculation^ of the force-constants. QE is a lin- 
ear response based method : the density functional per- 
turbation theory (DFPT).^ The dynamical matrix for 
the phonon excitation of a system is obtained from the 
ground state electron charge density and its linear re- 
sponse to a distortion of the ion-core geometry. We refer 
the reader to a recent article^ for further computational 
details on Fei-ajCr^; alloys. This alloy, being a basic in- 
gredient of stainless steel, is a technologically important 
structural material, dominated by force-constant disor- 
der and hence should serve as a critical test of our theory 
for the thermal transport properties. 

We find that the disorder induced scattering effects 
on the thermal conductivity, k(T), is relatively large in 
the low frequency regime. k(T) shows a quadratic T- 
dependence in the low temperature range, where only low 
energy vibrations are excited, and then smoothly rises 
to a T-independent saturated value at high T. Ther- 
mal diffusivity manifests the effect of disorder in a more 
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FIG. 1. (Color Online) (Top row) Dyson's Equation due to 
scattering diagrams for the single particle averaged Green's 
functions for disordered alloys. (Middle Row) Bethe-Salpeter 
equation for the response functions in disordered alloys. (Bot- 
tom Row) Key to diagrams. < G > is the averaged disorder 
renormalized Green's function, < k > is the two-particle cor- 
relation function related to Thermal conductivity, E is the self 
energy and Sotr is the disorder-renormalized effective current. 



dramatic fashion, and gives an idea about localization. 
Based on our calculation on Fei-ajCij; alloys, a large frac- 
tion (> 90 %) of vibrational eigenstates are found to be 
localized with the maximum localization near 50-50 com- 
position, where the disorder scattering is maximum, as 
expected. 

For disordered materials, the lattice thermal conduc- 
tivity requires the configuration average of the response 
functions of the kind (see Ref. 6), 

{{k{z,,Z2,T)))^ 

^Tr[((S(k,r)G(k,zi)S(k,r)G(k,Z2)))] , (1) 

where S is the heat current operator and G is the phonon 
propagator. (( )) denotes configuration averaging. 

The right hand side of Eq. (1) involves the configu- 
ration average of four random functions whose fluctua- 
tions are correlated. Unlike the conflguration average of 
a single particle Green function ((G(k, z))), which can 
be calculated via a perturbative self-consistent Dyson's 
equation (shown diagrammatically in the 1st row of Fig. 
1), the average of a two-particle correlation function such 
as that in Eq. (1) is non-trivial. The zeroth order ap- 
proximation for such an average is the one which as- 
sumes the fluctuations between all four random func- 
tions to be uncorrelated, and expresses the average of 
the product as the product of the averages (as in the 
so called Virtual Crystal Approximation (VGA)). The 
inherent correlation, however, requires the contributions 
from averages taken in pairs, triplets and all four random 
functions. Such disorder induced corrections can be cal- 
culated very efficiently within a Feynman diagrammatic 
technique (details have been discussed in Ref. 6), which 
renormalizes both the phonon propagators as well as the 
heat currents to provide a mathematical expression for 
((k)) with an effective heat current ScS related to the 
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FIG. 2. (Golor Online) Frequency dependence of Thermal 
conductivity and Joint density of states for Fe53Cr47 alloy at 
T — 200 K. Solid line shows the result including all disorder- 
induced corrections -I- the vertex correction (middle row of 
Fig. 1) and dashed line including the VCA average alone. 



self-energy of the propagators (shown by 1st diagram on 
RHS of the middle row of Fig. 1). The last term in the 
middle row gives the contribution from the so called ver- 
tex correction arising out of the correlated propagation. 
For a harmonic solid, thermal diffusivity has a similar ex- 
pression as ((k)) except the product of five random func- 
tions instead of four. A similar diagrammatic procedure 
has been used earlier by us® to calculate the configuration 
averaged thermal diffusivity as given by, 



, d^k 



Tr [((SmG(k,z/')S(k)3mG(k,i/')S(k)Q=mG(k,i/)))] . 

In Fig. 2, we display the frequency dependence of lat- 
tice thermal conductivity k(;/) and the scaled joint den- 
sity of states J(i') at T = 200 K for the Fe53Cr47 alloy. 
It is obvious from the figure that the transition rate r 
(related to the heat current operator) is strongly depen- 
dent both on the initial and final energies throughout 
the phonon spectrum i.e. k{v,T) ^ \t{v,T)\J(i'), where 
J{v) (shown by dot-dashed line in Fig. 2) is given by 



/ d^' j ^Tr[3?«((G(k,i.')))Sm((G(k,i.' + j/)))] . 

The effect of disorder-induced renormalized corrections 
(black solid lines) to the zeroth order virtual-crystal- 
approximation (VCA) (blue dashed line) thermal con- 
ductivity is not significant, and becomes negligibly small 
beyond v = 2.7 THz. The traditional single site mean- 
field approximation is, therefore, expected to describe 
well the multiple scattering phenomenon associated with 
the high frequency mode, deviating only in the low fre- 
quency range where the higher order corrections become 
important. Notably, both k.{i') and J{v) curve has a 
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FIG. 3. (Color Online) Temperature dependence of Thermal 
conductivity (k) for three Ve^-xC^x alloys. Inset shows the 
quadratic T-dependence of k in the low T-range. 



dip at a very small energy (i/ ~ 0). Such a dip reflects 
the missing intraband contribution (k^^) to the conduc- 
tivity. The origin of this dip is a natural outcome of 
a smooth convolution of two Green matrices G(k, i/') 
and G(k, i^' + v) (or two smooth DOS). A similar dip 
has also been reported by Feldman et al.'^ in amorphous 
Si and Sii_a;Gea; alloys. Unlike our case {k{v) — as 
1/ — > 0), this dip in their calculation stands at a finite 
value as u ^ Q. These authors have introduced an arbi- 
trary Lorenzian broadening of the delta functions in their 
Kubo-Grcenwood expression for k, while in our calcula- 
tion this arises naturally through the disorder induced 
broadening of the spectral function 3TO[G(k, v)]. An ex- 
trapolation of our K(i^)-curve (see Fig. 2) from a value 
just above = to a value at = yields an estimate 
of the dc thermal conductivity, which comes out to be 
24.7 W/m/K for the present Fe53Cr47 alloy at T = 200 K. 
Literature survey shows a lack of available experimental 
data for concentrated Fei-xCr^; alloy, however there ex- 
ist some on the dilute Cr-alloys.^" For example Kcxpt for 
X = 0.25% Cr is - 22 W/m/K, with which we shall com- 
pare our theoretical estimate below. 

Figure 3 shows the temperature dependence of ther- 
mal conductivity for three Fei-a-Cr^; alloys. Note that 
k(T) behaves quadratically (see inset) in the low tem- 
perature regime (T < 20 K) where only low-energy vi- 
brations are excited. As the temperature is increased fur- 
ther, the T-dependence of k becomes much milder and 
eventually reaches a T-independent saturated value. The 
origin of such a high T-saturation is not very well de- 
scribed by most previous theories. Within a harmonic 
approximation, such a saturation mainly arise from the 
T-dependence of Einstein specific heat piece of the con- 
ductivity expression.^ The intrinsic harmonic diffusion 
of higher energy delocalized vibrations are mostly re- 
sponsible for the relevant dominant mechanism in this 
regime. Another qualitative explanation can be that : 
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FIG. 4. (Color Online) Thermal diffusivity (top) and DOS 
(bottom) vs. phonon frequency {v) for Fe53Cr47 alloy. In- 
set shows an approximate linear i^-dependence of D[v) above 
2 THz. Vc locates the mobility edge and area under the shaded 
region gives an estimate of the fraction of delocalized states. 



the phonon-phonon scattering in this high T-range be- 
comes so strong that the phonon MFPs reach a minima, 
and further enhancing the disorder scattering by rais- 
ing temperature would not cause any further reduction 
in the MFP, hence resulting in a T-independent thermal 
conductivity. This, however, is just a physically plausible 
explanation based on the MFP and is not intended to re- 
flect a known outcome of the proposed theory itself. One 
can also notice an overall reduction of k with increasing 
disorder (x), as expected. Such effect usually reflect the 
scattering arising out of the difference in masses, radii 
and force constants between the host lattice atoms and 
impurities. In the present Fei-^^Cra, alloy, however, this 
scattering is mainly dominated by a large difference of 
force constants between Fe and Cr atoms in the alloy, 
while their masses and radii are almost similar. As far 
as the comparison goes, looking at the solid black curve 
of Fig. 3 for x=0.25% Cr, one can notice the saturated 
value of thermal conductivity (room temperature value) 
to be 22.8 W/m/K which is in good agreement with the 
experimental value of 22 W/m/K. 

Next, we examine the effect of disorder scattering on 
the vibrational eigenstates and hence the localization of 
the phonon modes based on a thermal diffusivity calcu- 
lation. In Fig. 4, we show the thermal diffusivity (top 
panel) and the phonon density of states (bottom panel) 
vs frequency for the Fe53Cr47 alloy. Above v ~ 2 THz, 
D(i^) decreases smoothly (approximately linear in v) with 
a critical frequency Vc = 3.55 THz, where D(i/) vanishes 
to within a very small level of noise. This regime is 
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FIG. 5. Mobility edge (top) and the percentage of localized 
states (bottom) vs Cr-concentration for Fei_a;Cra; alloy. 

shown, for clarity, as a log-log plot within the inset of the 
top panel. The calculated critical exponent a ~ 1.011 
agrees with the scaling and other theories of Anderson 
localization.^^ The critical frequency Vc locates the mo- 
bility edge above which the diffusivity strictly goes to 
zero in the infinite size limit, and the allowed vibrational 
states above (below) this frequency remain localized (de- 
localized). This is shown by the area under the shaded 
(unshaded) region in the bottom panel which gives an es- 
timate of the percentage of delocalized (localized) states. 

The location of mobility edge (j^c) and the percentage 
of localized states (calculated from the area under the 
curve in Fig. 4) with varying Cr-content in Fei-a^Cri: 
alloy is shown in Fig. 5. Such a non-monotonous vari- 
ation of the fraction of localized states is an artifact of 
the varying band- width of the phonon spectrum with x. 
Unlike the mass dominated Sii-j-Ge^; alloys^ which show 
an increasing percentage of localized states, towards the 
upper end of the phonon spectrum, with increasing Ge- 
concentration, the Fe-Cr alloys show maximum localiza- 
tion at xcr = 47 %. Such a maximal localization may 
be due to the dominance of the force constant disorder 
in the present alloy which causes an enhanced disorder 



scattering at a; = 47 % and hence localize the vibrational 
modes maximally. 

An alternative way of estimating the fraction of lo- 
calized (delocalized) states is to calculate the so called 
^'inverse participation ratid^ ^/Pj defined as, l/pj = 
J^f^Iid^k/S-K^) ei;(k), where e^(k) is the /xth Carte- 
sian component of the normalized polarization vector of 
the 7th mode, p-y measures the number of atoms on 
which 7th vibrational mode has significant amplitude. 
l/p-y — > for delocalized mode, but remains finite for lo- 
calized modes. Although this procedure provides a quick 
assessment of localization, it suffers from a shortcoming 
which arises quite often from the existence of an unex- 
pected few localized modes in the low frequency regime 
as discussed earlier. ^'^^ This is mainly due to a sensitive 
dependence of p^ on the boundary condition used in the 
concerned model. In other words, finite-size theory (even 
for large model systems) cause an unphysical gap at the 
bottom of the spectrum, and the states closest to this 
gap may be or appear to be localized. The same states 
in a macroscopic sample, however, may not be localized 
but propagating (or may be resonant) .'^'-'^^ The percent- 
age of localized (de-localized) states calculated using the 
area above (below) Vc of DOS curve and using p.^ may 
differ, if calculated from such finite-size theories. How- 
ever, being a k-space based formulation, our theory does 
not suffer from such differences and is free from the unex- 
pected errors arising from the existence of few localized 
modes in the low energy regime. 

In summary, we combine a generalized Kubo- 
Grecnwood type formula with the linear-response based 
QE calculation to make a first principles prediction of the 
thermal conductivity and diffusivity of disordered Fe-Cr 
alloys. The effect of disorder-induced scattering on k is 
found to decrease with increasing phonon energy. Ther- 
mal conductivity shows a quadratic T-dependence in the 
low T-regime, increasing smoothly to a T-independent 
saturated value at high T. Thermal diffusivity provides 
an estimate of the location of mobility edge, which sub- 
sequently gives an idea about the disorder-induced local- 
ization in the system. Vibrational modes in the present 
Fbi-xCtx alloy are maximally localized at a; = 47 %, 
where the effect of disorder scattering is maximum. 
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